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Abstract. In this paper, we show that the infinite generalised Stirling matrices associated with 
boson strings with one annihilation operator are projective limits of approximate substitutions, 
the latter being characterised by a finite set of algebraic equations. 



1. Introduction 

The series of papers [H [21 [3] had two sequels. First one, algebraic, was the construction of a 
Hopf algebra of Feynman-Bender diagrams |10i [TT] arising from the product formula applied to 
two exponentials. Second one was the construction and description of one parameter groups of 
infinite matrices ^i9j and their link with the combinatorics of so called Sheffer polynomials. The 
object of this paper is to continue the investigation of those one-parameter groups by highlighting 
the structure of the group of substitutions. 

It is shown here how we can see this group as a projective limit of what will be called approximate 
substitution groups. 

First, we consider Boson creation and annihilation operators with the commutation relation 



Recall that [T| 

• the annihilation operator a, in the second quantization, represents an operator which 
changes each state lA*") of the Fock space (containing > 1) particles to another containing 
— 1 particles. One has here 



E-mail: 



[o, af ] = 1 



(1) 



a\N) 




(2) 



• the hermitian conjugate of the annihilation operator is the creation operator which 
changes each state |A^) of the Fock space containing particles to another containing 
+ 1 particles. One has then 

a^\N) = ^/]vTT|iV + l) (3) 

Starting from the vacuum |0), we can reach all the normalized states in the Fock space. There 
are indeed given by: 

\N) = ^|0) (4) 

2. Normally ordered form 

The noncommutativity of annihilation and creation operators may cause problems in defining 
an operator function in quantum mechanics. To solve these problems, we have to find some 
suitable form which allows computing, reduction and equality test. One of the simplest and 
widely used procedure is the finding of the normally ordered form of the boson operators in 
which all stand to the left of all the factors a [l]. 

There are two well known procedures defined on the boson expressions: namely A^, the normal 
ordering and :: , the double dot operation, [ij. 

2.1. Normal ordering 

By normal ordering of a general expression F{a\ a), we mean AA[F(a1^, a)] which is obtained by 
moving all the annihilation operators a to the right, using the commutation relation ([T|). 

Example 2.1 Let w € {a,a^^}* a word given by w = aa^aaa^a. The normal ordering of w is 

aa'^ aaa'^a = (1 + a^a)a(l + a^a)a = + + aa^a^ + o^o^a^a^ 
= + + {I + a)a^ + a{l + a)a^ 



2a^ + 'ia'^ + a\l + a] 
2a^ + Sa^a^ + {a^fa'^ 



3 



Remark 2.2 Note that all coefficients of the normal ordering of a word (more precisely the 
coefficients of the decomposition of a word on the basis |(a^)-'a'|j are positive integers. This 
suggests that these integers count combinatorial objects fj^. 

2.2. Double dot operation 

The double dot operation : F{a'^,a) : is a similar procedure using another commutation relation 
i.e. [a, a^] = instead of = 1, i.e moving all annihilation operators a to the right as if 

they were commuting with the creation operators . 

Remark 2.3 The double dot operation :: is a linear operator which can be directly defined, for 
a word w E {a,a^}*, by: 

where \w\x stands for the number of occurences of a symbol x in the word w. 

Example 2.4 We take the same word as above: w = aa'^aaa'^a. 
The double dot operation gives 

: aa^aaa'^a : = a^a^aaaa 



3. Combinatorics of the normal ordering 

The BeU and StirUng numbers have a purely combinatorial origin [8], but in this communication, 
we will consider them as coefficients of the normal ordering problem jl] . 

The general word w G {a,a^}* with letters in {a, a^}, i.e. Boson string, can be described by 
two sequences of non negative integers r = (ri, r2, • • • , rjv/) and s = (si, S2, • " " j ■sm), so that we 
define 



(5) 



and d = (r^ — Sm), n = 1, • • • , M represents the excess (i. e. the difference between the 



m=l 



number of creations and the number of annihilations). Then, the normally ordered form of 
is given by 



(at)""'^5r,s(n,A:)(at)V, if d > 0; 
fc=0 

oo 

Sr,s{n, k){a^Ya'') (at)"l'^l, otherwise. 



k=0 



where quantities Sr,s{n, k) are generalizations of standard Stirling numbers OS]. 
The generalized Bell polynomials Bj-^sin, x) and the generalized Bell numbers -Br.sl'^) are defined 
respectively by 



5r,s(n, 1) 



^ Sr,s{n, k)x^ = 5r,s(n, k)x^ 

oo nr 



(6) 
(7) 



fe=o 



k=0 



Remark 3.1 Notice that 



1, for k = nr; 
0, for k > nr. 



Example 3.2 Using a computer algebra program, we get the following matrices. 



w = a^a , we get the usual matrix of Stirling numbers of second kind S{n, k)[S\ 



/ 1 

1 

oil 

13 1 

17 6 10 

1 15 25 10 1 

1 31 90 65 15 1 

V ; ; ; : ; ; 



w = a^aa^ , we have 
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18 
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1 











24 


96 


72 


16 


1 








120 


600 


600 


200 


25 


1 





720 


4320 


5400 


2400 


450 


36 


1 



w 



one gets 



1 
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4 


1 




















12 


60 


54 


14 


1 














144 


1296 


2232 


1296 


306 


30 


1 








2880 


40320 


109440 


105120 


45000 


9504 


1016 


52 


1 



Remark 3.3 In each case, the matrix {Sr,s[n-,k\)n>o,k>Q is of a staircase form and the dimen- 
sion of the step is the number of a's in the word w. Thus, all the matrices are row-finite and 
are unitriangular iff the number of its annihilation operators is exactly one. Moreover, the first 
column is (1, 0, • • • , 0, • • •) iff w ends with a (this means that J\f{uj^) has no constant term for 
alln> O; igU/. 

A word w with only one annihilation operator can be written in the unique form w = 
{a'<Y-Pa{a^)P. Then: 

• if P = ^! is the matrix of a unipotent substitution. 

• if p > 0, is the matrix of a unipotent substitution with prefunction J^. 
4. Approximate substitutions 

We define here the space of transformation matrices and its topology, and then we concentrate 
on the Riordan subgroup [7] {i.e transformations which are substitutions with prefactor func- 
tions) [11[9]. 

Let be the vector space of ah complex sequences, endowed with the Prechet product 
topology [6^, 12]. The algebra £(C^) of all continuous operators — > is the space of 
row-finite matrices with complex coefficients (a subspace of C^^^). Let M be a matrix of this 
space. For a sequence A = {an)n>o, the transformed sequence B = MA is given hy B = {bn)n>o 
with: 



bn = y^M(n,fc)afc. 

A;>0 



(8) 



We will also associate to A (see paragraph ()4.2|) ) its exponential generating function 



n>0 



4-1- Substitutions with prefunctions 

We now examine an important class of transformations: the substitutions with prefunctions. 
We consider, for a generating function /, the transformation 



g{x)f[(j){x) 



(10) 



where g{x) = 1 + X]5S=i 9nx'^ and (/)(x) = x + J2'^=2 4> 

arbitrary formal power series. The 
mapping 7^^^ is a linear application [9], the matrix of this transformation Mg^^p is given by the 



transforms of the monomials ^ hence 



n>0 



X 



(11) 



Remark 4.1 If f{x) = e^^ , Eq. comes down to the Sheffer condition on the matrix ofT. 
It amounts to the statement that 



V r(n,fe)-/ = 5(x)e^'^(^) 
^-^ n! 

n,fc>0 



(12) 



where T G £(C^) is a matrix with non-zero two first columns. 



4-2. Approximate substitutions 

In this section, we define approximate substitutions matrices and give a way to determine 
whether an unipotent (lower triangular with all diagonal elements equal to 1) matrix is a matrix 
of an approximate substitution. 

Definition 4.2 Let M € iC^o-n]xlo--n] ^ unipotent matrix, M = M[i, A:]o<i,fc<n = (flifc)o<j,fc<n/ 
is called matrix of approximate substitution if it satisfies the following condition: 



where 



Ck 



M 



Co 



A:! 



for all 



Ck = Y,M[i,k]- 

/ 1 

ai,o 1 

02,0 02,1 1 

03,0 03,1 • • • 1 

04,0 04^1 • • • a4^k • • • 

05,0 05^1 • • • Os^fc • • • 

V ; ; : : ; 



0<k<n 



(13) 



Co Cl Cfc 



Thus Ck represents the exponential generating series (here a polynomial) of the k^^ column 
(hence cq, ci are respectively the exponential generating series of the 1** and the 2"'^ column) 

and is the truncation, at order n, of a series. 

J n 

We consider now the set of matrices with complex coefficients noted by C^^^ [3] and 
let c[°" "'^[° ""'l be the set of all matrices of size (n + 1) x (n + 1). Let also r„ be the 
truncation of the matrices taking the upper left principal submatrix of dimension (n + 1), hence 

rn{M) = [M[i, k] ) . Thus, we get a linear mapping 

V / 0<i,k<n 



■n] X [0---n] 



It is clear that r„ is not a morphism for the (partially defined) multiplication (i.e. 
rniAB) / rniA)rniB) in general). 

We consider now £T(N, C) the algebra of lower triangular matrices and CT{[0 ■ ■ ■ n],C) the 
matrices of size ([0 • • • n] x [0 • • • n]) obtained by the truncation r„. Then 

Tn : £T(N,C) — > £T([0---n],C) 

and, this time, Tn preserves multiplication (r„ is a morphism). One has the diagram 

— ^ (p[0---ri]x[0---n] 



J2 



£T([0---n],( 



(j-<NxN : 

Ji 

£T(N,C) 

where and ^2 are two canonical injections. 
Remark 4.3 We can write 

£T{N,C) =lim {£T{[0 ■ ■ ■ n],C)) 
which means that £T([0- • •n],C) is the projective limit o/£T(N, C) 



(14) 



4.3. Random generation 

Our motivation, in this section, consists in approximating the matrices of infinite substitutions 
by finite matrices of (approximate) substitutions. We are then interested in the probabilistic 
study of these matrices. To this end, we randomly generate unipotent (unitriangular) matrices 
and we observe the number of occurrences of matrices of substitutions. 
The construction of unipotent matrices is done as follows: 

(i) Create an identity matrix. 

(ii) Fill randomly the lower part of the identity matrix (strictly under the diagonal) using 
numbers which follow a uniform law in [0, 1]. 

(iii) Multiply those numbers by the range previously chosen. 

(iv) Test if the built unipotent matrix satisfies Eg. (1131) . 

We start by giving some examples of our experiment which are summarized in the table below: 



bize 


IN umbers or drawing 


Range of variables 


Probability 


[3 X 3] 


300 


[1 • • • 10] 


1 






[1 • • • 100] 


1 






[1 • • • lOOOOJ 


1 


[4 X 4] 


275 


[1 ■ ■ ■ 10] 


0.0473 






[1 ■ ■ ■ 100] 


0.0001 






[1 • • • 10000] 


0+ 


[10 X 10] 


1500 


[l-.-lO] 


0.0327 






[1 • • • 100] 


0+ 






[1 • • • lOOOOJ 


0+ 



According to the results obtained, we observe that the substitutions matrices are not very 
frequent. However, in meeting certain conditions such as size, the number of drawings and the 
range of the variables, we can obtain positive probabilities that these matrices appear. 
Let us note that the smaller the size of the matrix the more probable one obtains a matrix of 

substitution in a reasonable number of drawings. 

We also notice that, if wc vary the range of variables, and this in an increasing way and by 
keeping unchanged the number of drawings and size, the probability tends to zero. We also 
notice that the unipotent matrices of size 3 are all matrices of approximate substitutions. This 
is easy to see because the exponential generating series of the 3'"'^ column will always have the 

form Cfe = 

Thus, we can say that the test actually starts from the matrices of size higher or equal to 4. 

Result 4.4 Let r represent the cardinality of the range of variables and nx n be the size of the 
matrix. 

According to the results obtained; we can say that the probability Pn of appearance of the matrices 
of substitutions depends on r and n and we have the following upper bound: 



Pn < (15) 
7" 2 

which shows that 



Pn — ^ as n — >■ oo (16) 

Conjecture 4.5 One can conjecture that the effect of the range selection vanishes when n tends 
to infinity. More precisely: 

1" 2 
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